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Abstract. New exact solutions emerge by replacing the dust source of the 
Lemaitre-Tolman-Bondi metrics with a viscous fluid satisfying the monatomic gas 
equation of state. The solutions have a consistent thermodynamical interpretation. 
The most general transport equation of Extended Irreversible Thermodynamics is 
' satisfied, with phenomenological coefficients bearing a close resemblance to those 

characterizing a non relativistic Maxwell-Bolzmann gas. 

\Q ■ 1. Introduction. 

(N . 

Dust is the prefered matter source in general relativistic models if rest mass is 
the dominant form of matter-energy! 1 ' 2 !. Exact solutions with a dust source might 
provide an appropriate description of the dynamics of astrophysical systems, but 
lack any information on their thermal properties and evolution. A more complete 
description can be achieved by the replacement of the dust source with a non- 
relativistic monatomic ideal gas. Since natural processes are irreversible and ideal 
gases (in general) are characterized by non-zero transport coefficients, a hydrodynamic 
model of this type of source should be based on the momentum energy tensor of an 
imperfect fluid. The purpose of this paper is to present a new class of exact solutions 
in which such matter tensor, subjected to the ideal gas equation of state, is the source 
of the Lemaitre-Tolman-Bondi (LTB) metrics' 1 ' 2 ! . However, Einstein field equations 
applied to these metrics (in comoving coordinates) necessarily restrict the fluid tensor 
to zero heat conduction and a non-accelerating 4-velocity, so that shear viscosity, 
associated with an anisotropic pressure tensor, becomes the only dissipative flux and 
the resulting shear viscosity tensor must be of a very special kind: its divergence 
must exactly ballance a non-zero gradient of the equilibrium pressure. Since, in 
general, classical ideal gases have non-zero heat conduction and pressure gradients 
produce non-zero 4-acceleration, the lack of these features obviously limits the physical 
applicability of the models, but is justified by the mathematical simplification of the 
field equations. Fortunately, though, the above mentioned limitations are compensated 
by the obtension of exact solutions that are thermodynamicaly consistent: the models 
derived in this paper are wholy compatible with the requirements (when shear viscosity 
is the only dissipative flux) of theories dealing with irreversible phenomena, such as 
Extended Irreversible Thermodynamics' 3-10 ' and Kinetic Theory' 3 ' 11 ' . Bearing in 
mind their shortcomings, these fluid models of a classical ideal gas, evolving along 
adiabatic (zero heat conduction) but irreversible (non-zero viscosity) thermodynamical 
processes' 11 ' 12 ' can be considered a reasonable initial improvement on dust sources in 
the description of astrophysical systems in classical conditions. 
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2. Monatomic ideal gas sources. 

A non-relativistic, non-degenerated, monatomic ideal gas is characterized by the 
equation of state' 7 ' 9 ! 

3 

p = mc 2 n+ -nkT, p = nkT (1) 

where p, n,p, T are matter-energy and particle number densities, equilibrium pressure 
and temperature, m, k are the particles' mass and Boltzmann's constant. As 
mentioned in the previous section, the hydrodynamical description of this ideal gas 
will be provided by solving Einstein's field equations for the momentum-energy tensor 
of a fluid with dissipative fluxes 

T ab = pu a u b +ph ab + U ab + 2c-y a ii b \ (2) 
h ab = a' 2 u a u b + g ab , u a lT b = 0, II a a = 0, u a q a = Q, 

where the equilibrium variables p,p satisfy (1), u a is the 4- velocity of the particle 
frame' 5 !, IT ob is the shear viscous pressure and q a is the heat conduction vector (bulk 
viscosity is negligible for a classical ideal gas' 3 ' 9 !). This matter tensor must satisfy 
energy-momentum ballance, as well as matter conservation and entropy production 
laws 

T ab :b = (3a) 
(nu») ;o = (36) 
(n S u a \ a > (3c) 

where s is the entropy per particle, which together with q a and 11^ must comply with 
suitable transport and constitutive equations of Irreversible Thermodynamics. 



3. A class of exact solutions. 

Consider the Lemaitre-Tolman-Bondi (LTB) metric ansatz associated with spherically 
symmetric dust solutions! 1 ' 2 ] 

ds 2 = -c 2 dt 2 + Y 'l dr 2 + Y 2 (d9 2 + sin 2 6dcj> 2 ) (4) 
1 k j 

where k = 0, ±1, Y = Y(t, r), / = f(r) and Y' = Y r . Following the strategy outlined 
previously, we shall consider the fluid tensor (2) (instead of dust) as the matter source 
of (4). If the spacial coordinates in (4) are comoving, we have u a = c6", and so, the 
4- velocity is a geodesic field (u a . b u b — 0) characterized by two kinematical invariants: 
the expansion scalar 9 = u a . a and the shear tensor a ab = M( a:b ) — (l/3)Qh ab , given 
explicitly by 



_ 2Y Y' 



(5a) 
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a\ = diag[0,-2<7,a,<7], a ab a ab = 6a 2 , a = ± ^ - ^ (56) 

where y = ?i a y a = cy t . The field equations for (2) and (4) imply vanishing heat 
coinduction: q a = q(t, r)5 r a = 0, while the shear viscosity tensor in (2) (a traceless 
symmetric tensor) takes the form 

n Q b = diag [o, -2P, p, p] , n afc n afc = 6P 2 (6) 

where P = P(t, r) is an arbitrary function to be determind by the field and transport 
equations. If q a = u a = the ballance equation (3a) yields 

p+(p + P )e + a ab Il ab = => p+\Qp + AaP = (7a) 



n 



h cd ) = o 



bc:d 



-p' + IP' + 6Py = 



(76) 



where p — u a p a , p = u a p a , the quantities a, P are defined in (5) and (6), and the 
relation n + n<d = from (3b), as well as (1), were used to derive the expressions in 
the rhs of the implication sign. 

Einstein's field equations for (2), (4) and (6), together with the integration of the 
matter conservation equation (3b), yield 



Y 2 = 



8ttG 

3c 2 r 



M + W 



-V 2 / 2 



where: M = mc 2 J n Q Y^dr, W =^j n o kT o Y S Y u dr 



(8) 



(9) 



and the subindex " " below any function (as in n ,T ,Y ) will denote henceforth 
evaluation of the function along a suitable intial hypcrsurface t = t . With the help 
of (8) and (9), the thermodynamical variables and a are given by 



yyi 

YY 
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The solutions characterized by (1)-(15) become fully determined once (8) is integrated 
for specific initial conditions prescribed by selecting n , T , Y 0l f (one of these functions 
can always be eliminated by an arbitrary rescaling of r). Ballance equations are 
identicaly satisfied, this follows by substitution of (5), (11), (12) into (7). Two 
important limiting cases emerge: (a) Lemaitre-Tolman-Bondi dust solutions are the 
zero temperature limit (T = 0^>T = p = P = 0), and (b) The shear- free particular 
case (a = 0, for all fluid worldlines) is the FLRW metric characterized by /'/ / = Y'/Y, 
Y = R(t)f, Y" = R f, rip = Tq = 0, whose source is a perfect fluid (P = 0) satisfying 
the equation of state (l)' 13 '. 



4. Thermodynamical consistency. 

After integrating the field equations (pending the integration of (8)), it is necessary to 
examine the thermodynamical consistency of the viscous fluid source. For such a fluid, 
Extended Ireversible Thermodynamics associates a generalized entropy current' 3-10 ! 
satisfying (3c) and relating the deviation from equilibrium due to viscosity 

(nsu a ). a > s > 0, s = s (e) - aTT ab n ab (16) 

where follows from the integration of the equilibrium Gibbs equation (Tds = 
d(p/n) +pd(l/n)) and a is a phenomenological coefficient. The evolution of the 
viscous pressure is in turn described by the transport equation' 8-10 ! 



rU cd h c a h d b +n ab 



2 '\T V 



where Ii ab = II ab . c w c , and the phenomenological quantities t, 77 are the relaxation 
time and the coefficient of shear viscosity. Equation (17) is the most general transport 
equation available when shear viscosity is the only dissipative effect. In order to deal 
with a mathematicaly simpler problem, less general equations are often suggested. 
For example, Eckart's equation (r = 0) or the "truncated form" without the term 
involving u c and T in (17). However, the former leads to acausal propagation equations 
and unstable solutions' 4,9 ', while numerical studies in FLRW spacetimes and other 
theoretical arguments' 8-10 ' point to unphysical effects in truncated forms of viscosity 
transport equations. 

In order to solve equations (16) and (17) for (1 )- (15) , the coefficients a, 77, r must be 
derived or specified. As a convenient reference providing a guideline on how to inferc 
physicaly reasonable forms of these quantities, consider a classical monatomic ideal 
gas characterized by a Maxwell-Boltzmann distribution near equilibrium' 3 ' . Grad's 
fourteen moments method provides the following forms for a, 77, r for this ideal gas 



k_ 

7o n v ~ 7 ' ' 1 mb *-- MB y-„ - MB 4r] MB nT V 



^«b = — («). V MB =PT MB (b), a MB = —M°^ = — 2 (c), (18) 



where 7 is a collision integral (eq (3.41) of [3]). A reasonable test of the physical 
viability of the solutions follows from verifying if the forms of n, T, p, Ii ab and 
a ab in (5), (6), (11)-(15) are compatible with equations (16) and (17) in which the 
phenomenological coefficients are given by (18). This is examined below. 
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Substitution of (11)-(13) into (16), with a given by (18c) and obtained from 
the integration of the equilibrium Gibbs equation for (1), yields 



s = s + k In 



rp \ 3/2 



Y'/Y 



$3/2 



3, fn 

~k . 



(19) 



where s (r) is an arbitrary initial value of s. Regarding (17), if we assume: i] = pr, 
with p given by (12), then this transport equation is satisfied for r given by: 
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(20) 



While the form of r\ can be justified as being formaly identical to the Maxwcll- 
Boltzmann form (18b) (a remarkable fact), the form of r in (20) is acceptable as 
long as this expression behaves as a relaxation parameter for the system, namely: it 
should be a positive function and its relation with s (computed from (19)) must be 
such that it garantees s > along u a for all fluid worldlines. Idealy, of course, r should 
be related to the mean collision time of the gas and should also, somehow, relate to or 
approach the Maxwell-Boltzmann form (18a). Evaluating s from (19) and comparing 
with (20), we find the following simple and elegant relation between s and r 



HI) \ 



■u i — y - 

p J T 



(21) 



which implies: s > r > 0. This relation is in fact a direct consequence of the 
general expression: s = H ab H ab / (2rinT) (see [4], eqn. (18)), together with equations 
(3b), (3c) and (18b). Therefore, from equations (11) and (19)-(21), the following are 
necessary and sufficient conditions forT>0,r>0,s>0 



* > 

o-n < o 



(22 a) 
(226) 



while the conditions insuring that s decreases for increasing r (f > s < 0) follow 
from (20)-(21) 



f > 0, 
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(22c) 



The set (22) provides the necessary and sufficient conditions for a theoreticaly 
consistent thermodynamical description of the solutions, within the framework of 
equations (16)-(21). However, in order to work out these conditions in a more explicit 
and direct manner (especialy (22c)), we need to know the explicit form of Y'/Y from 
integrating (8). 
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5. The case k. 



o 



0. 



Consider the integration of (8). For the simplest case, k — (the cases k — ±1 are 
examined in [14]), this yields 




, , X3/4 
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<5 = W/M, y = y/Yjj and F is the elliptic integral of the first kind with modulus l/\/2 
and argument <p given by: cos ip = (y — \/S) /(y + VS). With the help of (23b) and the 
definitions of ^P, f2, a, it is possible to provide specific examples showing that conditions 
(22) are satisfied for physicaly reasonable initial value functions n ,T ,Y Q . Consider 
the following trial functions, analogous to empirical models of mass distribution in 
spherical gallaxies and clusters' 15 ' 16 ! 



1 + 



1 + (Yo/RoY 



T a =T 



1 + 



1 + (Y /R ) 2 



R r 



(24) 



where n ,T ,a,b are positive empirical constants, while R is a characteristic length. 
Let us assume that (24) provides the initial particle number density and temperature 



profile at the time when matter and radiation decouple, t a 



'dec 



of a spherical 



inhomogeneity of gallactic mass (10 gm) in an expanding FLRW background 
characterized by an intial profile: n ss 10 5 cm~ 3 and T Q w 3 x 10 3 K. Taking m to 
be a protonic mass and considering a red shift factor z « 10 3 (so that R « 10 20 cm), 
reasonable estimates of a, b in (24) are a « 10 3 and b « 10. As illustrated by figures 
(la), (lb) and (lc) below, we have positive and monotonously decreasing T and s, 
and monotonously increasing t, so that the evolution of a configuration of this type 
is consistent with the thermodynamical requirements (16)-(22). In fact, it is possible 
to show that conditions (22) hold for all trial functions of the form (24), provided a, b 
comply with the constraint 



2 L 1 

3 { 1+ b 



<l + i 
a 



(25) 



together with kT /mc 2 < 1, Y > Y , conditions satisfied by a large choice of empirical 
constants. For n Q ,T given by (24) and (25), we also have' 14 ': r cx 1/n for Y/Y 3> 1, 
a form of r qualitatively analogous to (18a) and to relaxation times of globular 
clusters^ 17 ]. Further extension and development of the work presented in this paper is 
being undertaken.! 14 ]. 
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Figure Caption. The figures display: s (Figure (la)), T (Figure (lb)) and log 10 (r) 
(Figure (lc)), ploted vs. arctan(r) and arctan(y) (for y > 1, where y = Y/Y ), for the 
spherical inhomogeneity described in the text. Notice that s and r are everywhere 
positive and they respectively decrease/increase along increasing y (roughly equivalent 
to the fluid worldlines). This example is not meant to be a model of a "real" gallaxy, 
but simply to illustrate that conditions (22) hold for physicaly reasonable initial value 
functions n ,T . These plots were done with the symbolic package MapleV. 
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